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Abstract 

We discuss several topics related to the notion of strong hyperbolicity 
which are of interest in general relativity. After introducing the concept 
and showing its relevance we provide some covariant definitions of strong 
hyperbolicity. We then prove that is a system is strongly hyperbolic with 
respect to a given hypersurface, then it is also strongly hyperbolic with 
respect to any near by one. 

We then study for how much these hypersurfaces can be deformed and 
discuss then causality, namely what is the maximal propagation speed in 
any given direction. In contrast with the symmetric hyperbolic case, for 
which the proof of causality is geometrical and direct, relaying in energy 
estimates, the proof for general strongly hyperbolic systems is indirect for 
it is based in Holmgren's theorem. 

To show that the concept is needed in the area of general relativity we 
discuss two results for which the theory of symmetric hyperbolic systems 
shows to be insufficient. The first deals with the hyperbolicity analysis of 
systems which are second order in space derivatives, they include certain 
versions of the ADM and the BSSN families of equations. This analysis is 
considerably simplified by introducing pseudo-differential first order evo- 
lution equations. Well posedness for some members of the latter family 
systems is established by showing they satisfy the strong hyperbolicity 
property. Furthermore it is shown that many other systems of such fam- 
ilies are only weakly hyperbolic, implying they should not be used for 
numerical modeling. 

The second result deals with systems having constraints. The ques- 
tion posed is which hyperbolicity properties, if any, are inherited from 
the original evolution system by the subsidiary system satisfied by the 
constraint quantities. The answer is that, subject to some condition on 
the constraints, if the evolution system is strongly hyperbolic then the 
subsidiary system is also strongly hyperbolic and the causality properties 
of both are identical. 
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1 Introduction 



Consider a linear first order evolution system of differential equations in _/?""*" ^ 
on a set of tensor fields , 

where Greek indices denote a set of abstract indices of tensorial nature, we 
denote the dimension of that space by m, and A^°'p, and B" 13 are fields on 

The textbook definition of strongly hyperbolic systems is f|16 j .|10 | .and |17|'): 

Definition 1 The above system is strongly hyperbolic if for any co-vector 
LOi in _R", the matrix A^" ptOi has purely real eigenvalues and is diagonalizable. 

For constant coefficient systems this is a necessary and sufficient condition 
for the system to be well possed in the sense. Indeed, assuming constant 
coefficients and taking the Fourier transform in i?" of the fields 0", which we 
shall call 0" the solutions of the resulting system of linear ordinary differential 
equations are: 

rit,Lu)^{e^'^'^^+^^Yl3^^iO,^) 

and so, from Planchcrcl identity it follows that the solution will be bounded 
in with a bound that depends only on the norm of the initial data if and 
only if 

is bounded independently of u. This is the case if and only if the matrices 
A'uji have purely real eigenvalues and are diagonalizable. 

A matrix A^p has all its eigenvalues real and is diagonalizable if and only 
if there exists a positive bilinear form Ha/3 such that Ha^A'^ p is symmetric. In 
our case, we then have that for each given uji there would exists such H{uji)ai3 
which would symmetrize A^^pcui, namely H {u!i)[a\i3\A^^ j]UJi = 0. 

If one can find a common symmetrizer, Ha/s, for all A^"' puji then we say that 
the system is symmetric hyperbolic. 

In order to understand the difference between these two concepts, note that 
any linear combination of symmetric bilinear forms remains symmetric, but 
a linear combination of diagonalizable (or symmetrizable) matrices does not 
necessarily remains diagonalizable (or symmetrizable), for instance, A + AB, 
where 

-^{11) -^{11) 

is diagonalizable for all A 7^ —1, but not for A = —1. This implies that for strong 
hyperbolicity it is not enough to check that the matrices corresponding to some 
basis of ijJi's are diagonalizable, one has to check this for every co- vector uji. 
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Most systems in physics arc symmetric hyperbolic. ActuaUy it is not so easy 
to step into a genuinely strongly hyperbolic system, that is one which is not 
also symmetric hyperbolic. In particular, if n = 1 all systems are symmetric 
hyperbolic. So here is an example: 

Exeimple 1 

dtEi = d^Wji-il + a)d^Wij+adiW/ 
dtWij = diEj+p/2eijd^Ej, (1) 

where Latin indices reffer to tensors on R^, eij denotes the Euclidean metric 
there. Here, = {Ei,Wij). This system is strongly hyperbolic for all values 
of the paremeters a and (3 such that afi > 0, but there is no symmetrizer 
independent ofuji for a > or (3 > This system, subject to some constraints 
and gauge condition, is equivalent to Maxwell's equations, in vector potential 
form with Wij := diAj. This is the analog of what it is done in general relativity 
to obtain first order systems out of Einstein's equations. ^ 

The above definition seems to depend on a particular choice of time coordi- 
nate or foliation on the space-time. If that were the case this would be a way of 
choosing a preferred time function among all possible, and we know, on physical 
grounds, this can not be the case. So we must find out a covariant definition 
and explore the freedom left in choosing time functions. This we do in the next 
section. 

The above definition of strong hyperbolicity can be extended to general 
quasi-linear systems, where now the fields A^"js, and B"^ depend also on the 
unknowns, (j)" . In this more general setting strong hyperbolicity is necessary 
for well posedness, but for sufficiency one also needs smoothness of the matrix 
which diagonalizes the system, H{(p'^ ,cOi)ai3 which in general depends on Wj, 
the space-time point and the imknowns. We prefer here not to include, in 
the definition of strong hyperbolicity, this smoothness condition for we do not 
know whether it is truly needed, and so fundamental or it is just a temporary 
technical requirement for the proof -in fact we now we do not need it in the 
case of constant coefficients. 

How do we extend this concept to the case of manifolds? If we assume 
the manifold to be para-compact then we can assume we can work on every 
coordinate patch, extend the functions and equations on them to i?""*"^ and 
stablish well posedness on each one of them using pseudo-differential calculus. 
If in doing so we can get uniqueness and finite propagation speeds then we can 
patch the solutions found on each chart to get a global solution. The use of 
Fourier transforms seems to indicate one is dealing with very specific properties 
which might be valid only for particular manifolds (where we can apply Fourier 
series or transforms). This is not the case, what one actually is looking is 
at the high frequencies limit of the solutions of the given evolution equation 

^ There are other approaches where also the Bianchi identities are used to produce such 
systems, but now the variables include the Weyl Tensor. 
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and checking that in that hmit they do not grow unbounded with a growth 
rate which increases with the frequency. The behavior of bounded frequency 
solutions is not relevant for well posedness. This intuitive reasoning is backed 
by actual theorems arising from the pseudo- differential calculus technique, but 
in this article we shall try to stay away from them and only mention the relevant 
results. In the third section we shall look at the above problem, namely that of 
the causality properties of strongly hyperbolic systems, to check whether it is 
possible or not to use this definition in non-trivial manifolds. 

The last two sections shall present applications of this theory to two problems 
in general relativity. In the first we show how to analyze the hypcrbolicity 
properties of systems which are second order in space derivatives, they include 
certain versions of the ADM and the BSSN families of equations. This analysis 
is considerably simplified by introducing pseudo- differential first order evolution 
equations. We show how to obtain first order pseudo-differential equations for 
these systems and then how to check for strong hyperbolicity on them. These 
results have already been presented elsewhere and here are only sketched as 
illustrations of the theory. Well posedness for some members of the latter family 
systems is established by showing they satisfy the strong hyperbolicity property. 
Furthermore it is shown that many other systems of such families are only weakly 
hyperbolic, implying they should not be used for numerical modeling. 

The second is the study of constrained systems, on those systems, besides 
the evolution equations, there are extra equations which constitute differential 
relations that the initial data must satisfy in order to produce a valid solution. 
For that to happens some integrability conditions must be satisfied, in particular 
some providing the time derivatives of the quantities expressing the violation 
of these constraints in terms of themselves and their first derivatives, so that if 
they are zero at a given time they remain so for a finite time. Thus, to stablish 
constraint preservation, an understanding of the hypcrbolicity and causality 
properties of the evolution equations of these constraint quantities is needed. 
This is also of relevance for understanding the initial-boundary- value problem in 
the presence of constraints. In the last section wc shall sec which properties are 
inherited on the subsidiary system from the original evolution system. It turns 
out that the solution here is in terms of strongly hyperbolic systems: Strongly 
hyperbolic systems have strongly hyperbolic subsidiary systems. Furthermore 
some causality properties also are inherited. There are examples of systems 
whose evolution is symmetric hyperbolic but whose subsidiary system is only 
strongly hyperbolic, showing that also for this issue the concept of strongly 
hyperbolicity is needed. 

2 Strong Hyperbolicity 

We want now to introduce new definitions that reduce to the above one when 
a foliation is chosen, but which do not depend on any foliation, but rather in a 
property of an open set of vectors. We shall prove latter that these definitions 
are equivalent. 
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Let 



be a first order system on a fiber bundle wliose fiber has tangent space of 
dimension m and whose base manifold is a space-time, Af , of dimension n+ 1. ^ 

Definition 2 We shall say that this system is A- strongly-hyperbolic if there 
exists Ua e T+Af such that: 

i) A"'" is invertible, and 

ii) for each (non-contractible loop on the corresponding projective space) k(A)q = 
Aria + bJa where tOa is any other vector not proportional to ria, 

dim{span{U xi^rK ernlA""" l3Ka{X)}}) = m. 

Remarks: 

1. Choose a vector such that f^Ua ~ 1 and define ha^ — Sj" — riat^ then 
the above equation becomes, 

A°"/3na<''(V0)f - -A^^'pha'iVcbf, + J". 

Since A"" is invertible we recuperate the textbook definition by mul- 
tiplying by its inverse. 

The diagonalizability of {A""" pUa)^^ A'''^ ^uji, follows from the definition 
above, for it implies that the above matrix has a complete set of eigenvec- 
tors. Furthermore all its eigenvalues are real. 

2. Again it is not true that if B is a diagonalizable and invertible matrix and 
A is diagonalizable, then B~^A is diagonalizable. The matrices of the 
example above are a counterexample. 

3. In examples (e.g. the first order version of the wave equation in more than 
two space-time dimensions) it is easy to see that not every Ua for which 
A'^'^jsria is invertible is a good candidate for property ii), in fact it seems 
that only a connected subset of them in P" (see item 5 below) comprises 
the whole set of possible ria's. 

4. This definition still depends on a particular Ua and we must prove that 
the set of allowed Ua's is open. Otherwise it is a very restrictive definition, 
not better than definition ^ 

5. We should think on the loops above on the projective space P", where 
any two vectors are identified if they are parallel. This is a compact non- 
simply connected space. The Ka{X) loops are not contractible to zero, are 
compact, and continuous with respect to X,na and LUa- This implies, 

^Here we are using the notation of references |^. 
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• The roots of the determinant of A'^ continuous functions of the 
arguments. So the set of points in P" where the kernel of (Aria + 
Wa)^" is non-trivial is a continuous manifold of dimension n — 1. 
Furthermore, since (A""rjna) -fLOb is diagonalizablc for any aj},, 
and perturbations of eigenvalues due to diagonalizable matrices are 
Lipschitz, such manifold is Lipschitz. This sub-manifold can have 
some self-intersections where two or more eigenvalues coincide. 

• The curves K{X)a ■= Xria +C0a pinches off transversally this manifold. 
This follows easily from the fact that if we move continuously Wq the 
roots Xi must also move continuously and can not disappear, (become 
complex), because they are a complete set and are all real, so they 
must just move along the real axis. This implies that if we slightly 
perturb Ua the new curve must also transverse all these manifolds 
without loosing any root. The problem is then to assert that one 
does not looses any subspace in the kernel. 

• If we move continuously ria the (real) zeroes can not dissappear (be- 
come complex) so at most can move towards infinity, but those are 
precisely the points where Ua is no longer invertible. 

• f{na,cOa) := rfzm(UAGi?i^ern{/f""^Ka (A)}) is a upper-semi-continuous 

function of its arguments. That is, given any point (?!„, Ub) there is a 
always a neighborhood V such that /{ria, uia) < f{na,^b) V {na,COa) S 
V 

We are now ready to prove the main result of this section: 

Proposition 1 If for some Ua properties (i, ii) hold, then they also hold for 
any other co-vector in a sufficiently small open neighborhood of Ua- 

Proof 1 It is useful to think the problem (and use the corresponding topology) 

in real projective space P", where two vectors are identified if they are parallel. 
In that case the image of the curves Xua + uJa d^e compact non-contractible 
curves and so we can define a topology between them in a trivial way. 

We assume properties (i, ii) hold for som,e given Ua, and -for contradiction- 
that they do not hold for some fia = na + s6a, e > with e sufficiently small. 
That means that at som,e root Xi, the matrix {Xifia + i^v)^",- has a kernel with 
lower dimensionality than the one at (Ajna+cJa)A° where Xi is the corresponding 
root for this curve. Consider now the curve given by Xua + [uJa + A^eJa). At 
A = A; this curve intersects the above one and so its kernel at this point has also 
lost dimensions compared with the point (Xina+i^a) Since properties (i, ii) hold 
for this new curve there must be at least another point, XjUa + {uja + Xie5a) at 
which the kernel has increased its dimensionality compared with the nearby point 
XjUa +i^a where the fi,rst curve has a non-trivial kernel. But this is impossible, 
since the dimension of the kernel is upper semi- continuous. 

Proposition 2 Let Ua a co-vector such that properties (i, ii) hold, let S{na) be 
the connected component of vectors satisfying (i, ii), then the boundary of S{na) 
consists of vectors rua where A"ma ceases to be invertible. 
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Proof 2 Consider a one parameter family of vectors Va{X) '■= ria + XSa- Let 
Ao be the least upper bound of {X\given any co-vector ujathe curve Xva{X) + oJa 
has m real roots (counting multiplicities)}. For X G [0, Aq) the roots are real, 
and so by continuity the roots of the curve at Xq must be also real, and so the 
only possibility to loose a root is that it goes to infinity, that means the Va{Xo)A'^ 
has non-trivial kernel. We need to see now that this is a boundary point, that 
is, that at no point before Aq the dimension of the direct sum of the kernels is 
not smaller than m. Assume, for contradiction this happens, namely there is 
< Ai < Ao being the least upper bound of the set of X 's for which the dimension 
of the direct sum of kernels is m. Since the Kernel is a upper semi- continuous 
function at the point where it lowers its value must retain its upper value. So 
at Va{Xi) the dimension of the direct sum is m, but in arbitrarily near points 
that sum is smaller. But this contradicts that the set for which that sum is m 
is open. Unless at Ai the matrix ceases to be invertible. 

We know (see for instance 10 ) that ^ is a diagonalizable matrix with real 
eigenvalues if and only if the following inequality holds: 

\{sI + A)-^\<K/\Im{s)\ VseC, Im{s) ^ 

Thus, if a system is strongly hyperbolic for Ua, then, since A'^Ua is invertible 
and {A'^na)~^A'^uJa is diagonalizable with real eigenvalues. Thus we have, 

\{sA''na + A-uJa)-'\ = \{sI+{A-na)-^A-u;arHA''na)-^\ 

< Wn^r'WisI + {A''n,)-'A-u;a)-'\ 

< \{A-na)-'\K/\Im{s)\. (2) 

The reverse inequality is also true, thus, the above inequality is true if and 
only if A := (A°-na)~^ A°-u!a is diagonalizable and all its eigenvalues are real. So 
we can use also this property as an alternative definition of strong hyperbolicity. 

Definition 3 The above system is B-strongly hyperbolic if there exists Ua 
such that for any uJi and any s & C with '^s we have, 

lisA'^na-^A'^uJar^l < K/\Imis)\. (3) 

Proposition 3 A-strong hyperbolicity and B-strong hyperbolicity are equiva- 
lent. 

Let us see that this is indeed the inequality needed for proving well posedness. 
We consider the constant coefficient case and take a constant flat foliation, make 
Fourier transform in space 

^{t,uj,)^-^ J e-'^^^'<l>{t,x^)dV 
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and then a Laplace transform in time, 

/•oo 

4){s,uji) = / e^'''^4>{t,u!i)dt JRs > a > 0. 



Thus Eq. becomes, 
Thus we have, 

^(s.LOi) = (A"s + ^^■'■cc;,)-lA°^(0,c^^) 
Using Plancherel identity we get, with rj > a, 

T 2r,T poo 

< I \\{A"{rJ + ^C) + ^A^cu,)-'A"4>{0,u;l)\\'d<:du;- 



i^"; ' Joo 

g2jjT /-oo 



< ^^J^ \\A'4>{0,u;MA%rJ/u; + ^O+^A^Lu,)-'\\'LO-'dCdLo- 

^2 277T fOC 
2r,T 

where in the last inequality we have used inequality namely, 
\\{A°{r]/uj + iC)+iA^u;j)-^\\^ <Kuj/r] V UJ^ e S"-\ 

An alternative proof of proposition ^ is the following: assume inequality (jSJ 
holds for some Ua and all uja and try to probe it for Ha = ria + s6a where e is 
taken small enough. 

We have: 



A^'isha+uJa) = A''{sB.{na+e5a)+uJa)+isiA''{na+eSa) 
= A"n„(Trr-i + isiil + eiA^'nar^A'^Sa)) 
= A''naT{T + isi{I + eT-\A''na)-^A''6aT))T-^ (5) 

Here TTT~^ = {A^nay^isRl + ^"(wq + e6a)) where T is a diagonal matrix 
with all entries real. Such decomposition is possible because we have applied the 
assumption that {A"-na)~^A'^LOa is diagonalizable and has purely real eigenvalues 
for uja = t-^a + sSa. T and T depend on sj^, but since T only depends on the 
direction of sniua + e5a) + and the set of possible directions is a compact 
set, then \T\ and |T~^| are bounded. 
Thus, 

{A'^{sha+uJa))-^ = T{T + isi{I + eT-\A'^na)-^A'^5aT))-^T-\A''na)-^ 
= —T{y{isil + T)-\isiY{-irB^-^}T-\A%ia)^) 
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where B = eT~^{A"'na)~^A"'5aT But since F is diagonal and real we have that 
|s/(zs// + r)^^| is bounded and so for sufhciently small e the series converges 
and so the norm of the resolvent is bounded by some constant times l/|/m(s)|. 

The second proposition can also be proven using Eq. 

Consider the following alternative definition (due to R. Geroch T): 

Definition 4 A system is C-strongly hyperbolic if: 

i) For each Ua there exists hi3^{L0a) such that Apa ■= hp^{LUa)A°'^ a^^a is 
symmetric, and 

a) For some ria, and any uji,, h/3j{LJi,)A'^^ is symmetric and positive 
definite. 

We want to show that this definition is equivalent to the one already intro- 
duced. 

Proposition 4 A-strong hyperbolicity and C-strong hyperbolicity are equiva- 
lent. 

Proof 3 Assume the system is A-strongly hyperbolic. The kernel condition im- 
plies that for any LOc (A'^Ua)^^^ ^A'^'^ has a complete set of eigenvalues, thus 
there exists a positive bilinear form HspitOc) such that Hsp{ijJc){A°'na)~^^ ^A^^ a^i, 
is symmetric. But then, taking hsa{uJc) '■= F[sfj{uJc){A°'na)~^^ a we have 
hsa{L^c){A''°'jnb = HsyiiVc), and hsj{LOc)A'''^ a^^b = Hsfj{ujc){A°'nay^'^'yA'''^ a^^t 
and so both are symmetric and the first is positive definite. 

Assume now we have a C-strongly hyperbolic system. We want to show that 
for each (non-contractible loop on the corresponding projective space) K{X)a = 
Xua + LOa whcre to a is any other vector not proportional to Ua, 
dim{span{U\£iiKern{A'^°' pKa{X)}}) — m 

Now, u" is in such kernel for some A if and only if XA""" pUaU^ = A""" puJaU^ . 
Multiplying both members by h(ujc)-ya we get, 

Ajsu^u^ = XHysu^u^ 

with both Ay a, and H^p symmetric and the latter positive. But then the eigenvalue- 
eigenvector pairs are obtained as minimizing vectors of the function F{u°') := 
Ajgu'^ / Hjgu'^ on a diminishing sequence of H orthogonal subspaces, from 
which we conclude they form a complete set of vectors. 

3 Causality 

An important property of evolution systems is the propagation speed of high 
frequency perturbations. In particular the highest speed in any given direc- 
tion, via Holmgren's Theorem (see for instance ^7]) determines the domain of 
dependence of the system for any given initial hypersurface. 

■^Strictly speaking this is only asserted for linear systems with analytic coefficients, but will 
be generalized here. 
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Given a small enough region with local coordinate system {t, x') we consider 
fields of the form, = u^e"''^*""'^ , where is an eigenvector (and cr^^ the 
corresponding eigenvalue) of A^" puji. Then in the limit of high frequencies these 
fields are very close to solutions of the evolution equations and have velocities 
given by = Jj-. The surfaces tangent to a^dt—uji are called the characteristic 
surfaces of the system and along them the respective eigen-modes propagate, 
being the integral curves of V"" — (9f )" + v°- . The highest of these velocities for 
any given (which corresponds to the highest eigenvalue, since the dependence 
of the eigenvalues with respect to uji is homogeneous of degree one) determines 
the causality properties of the system. 

The argument for this last assertion is as follows: Consider a quasi-linear 
strongly hyperbolic system, say (|ll|l . and take a initial slice S which can be ei- 
ther continued to i?" or T", so that we can define on it either Fourier transforms 
or Fourier series and so introduce the pseudo-differential calculus needed to show 
well posedness for strongly hyperbolic systems. Extend the whole space-time 
to include a finite neighborhood of that extended slice, and give in the ex- 
tended slice initial data, (p^ with enough smoothness and decay properties as 
needed to use the existence theorems for strong hyperbolic systems (assuming 
the symmetrizer is smooth in all variables) to show that there exists a finite 
neighborhood of that initial surface where a unique solution, 0", exists. Con- 
sider now linear perturbations of that solution. Quasi-linearity implies they will 
satisfy an equation of the form, 

L"^<S0/3 at<5r - A(cj>, t, xy^pdM^ + B{(t>, t, xrpScj)". (7) 

That is, with the same principal part as the original one. 

We want to show that the domain of dependence of these perturbations is 
characterized by the highest eigenvalues for each direction. 

Following the same strategy as for Holmgren's theorem ^7j, we take a point 
p G S and two global coordinate systems with origin at p, (i,x') and {i,x^) 
which differ on the time coordinate in a cuadratic way near p so that we have 
two time foliations with the intersection of the past of one from one fixed time 
on with the future of the other form another fixed time on form a lens shaped 
region = > 0}n{£ < T}. We shall call the past bounary Eg and the future 
one Er- Choose the perturbation initial data, S(j)Q to vanish in Eq, but let it 
be otherwise arbitrary (but sufficiently smooth) everywhere else. The domain 
of dependence of such region will then be the region of space-time where the 
perturbation remains null, regardless of the initial values S4>q takes outside that 
region. If T is choosen sufficiently small then one can arrange for the normal to 
the t constant slices, fia = di, to be very close to the normals to the t constant 
slices, Ua — dt, so that with respect to both normals the evolution systems is 
strongly hyperbolic. 

We now consider solutions of the adjoint of the operator to lO, L"^ to 
the past of the surface i —T with vanishing initial data but with sources, 

K% = p\ 
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where p'a are a set of basis vectors times a function which inside the region fix 
is a polynomial of degree i in the coordinates and outside it decays exponentially 
to zero in such a way as keeping these sources smooth. Since the adjoint operator 
is also strongly hyperbolic the proof of well posedness for these systems will hold 
for each one of these source functions, so we have a dense set of solutions {V'^} 
satisfying: 

where the boundary terms go away since cji^ vanishes at the past boundary Eq, 
and V'7 vanishes at the future one, namely St- Thus, since the set of polynomials 
are dense, in we conclude that 0^ = inside the foliation. 

The result can be extended to the full solutions of the quasi-linear system 
by proving it at each step of the iteration: 

= A(</)„,i,x)'"^a,^^+i +i?(^„,i,a;)"05</)^+i. (8) 

Along the iteration the region where the solutions to the corresponding adjoint 
equations are valid can be shown to contain a given finite time interval which 
depends only on the initial data for 0". 

Since the lens shaped regions can be extended until the normal makes 
A"'°'f3ha singular we have the following result: 

Proposition 5 Given a first order quasilinear strongly hyperbolic system of 
equations the domain of dependence of regions is determined by the highest char- 
acteristics speeds along each direction. 

A different argument leading to causality can be obtained for constant coef- 
ficient systems by directly constructing their solutions as limits of grid solutions 
to their discretized versions. Indeed, as argued by Kreiss 9 , after embedding 
the system in or one can discretize the system using the method 

of lines and a centered difference operator of some given order. Such operators 
involve a finite number of grid points. ^ Thus obtaining a big system of lin- 
ear ordinary differential equations. It is easy to check the eigenvectors form a 
complete set and that the eigenvalues have no positive real part. Thus, for suf- 
ficiently small time-step size (depending on the space inter grid size) the system 
is stable and a solution exists and is bounded for any finite time interval. But 
at given grid size, the value of the solution at a given position in space and time 
depends, by construction, only on the values of the initial data grid functions in 
a finite region of the initial surface. When increasing the resolution the number 
of grid points which enter in such calculation increase, but not the size of the 
region of dependence, thus in the limit we see causality holds. 

Thus we see that strongly hyperbolic systems share many properties with 
the subclass of symmetric hyperbolic ones: 1.- One can prescribe initial data 
on any Cauchy surface, that is a surface such that its normal Ua satisfies the 

*We refer to a for the basic results we are quoting. 



11 



conditions on any of the equivalent definitions of strong hyperbolicity we have 
introduced; 2.- The causahty properties and propagation speeds of high fre- 
quency perturbations are derived from the characteristic surfaces of the system, 
namely surfaces at which the normal is such that A°'^ana is not invertible. 

But, are really necessary strongly hyperbolic systems in physics? It seems 
that all classical evolution equations in physics can be described in one way or 
another by symmetric hyperbolic systems, so why complicate matters introduc- 
ing strongly hyperbolic systems? In what follows we shall consider two examples 
where introducing strongly hyperbolic systems lead to new insights in physics. 

4 Well posedness of systems which are second 
order in space derivatives 

Here we sketch the analysis of well posedness for systems which are first order 
in time derivatives and second order in space derivatives. Details can be found 
in Main examples of such systems are the ADM equations of general 

relativity, |2] and the more recent modification of it called BSSN, fT5j . 

Since well posedness is a property related to the behavior of high frequency 
solutions in establishing it one looks at the equations and their solutions in 
that limit. It is a property of well posed systems that zero order terms do not 
change the properties of well posedness in the sense that if a system is well 
posed, then under the addition of a zeroth order term the system remains well 
posed 20). Thus, is we have a first order system we know that all we have to 
do is to look at its principal part, namely the terms containing first derivatives. 
So the strategy here to analyze these second order systems is to transform them 
into equivalent first order systems and then look at algebraic properties of their 
principal part. There are several ways of obtaining a first order system out of 
these second order ones. One of them, already used in [5], and is to add as 
variables all first order derivatives and look at the resulting larger system, where 
second derivatives are substituted by first derivatives of the added variables and 
new evolution equations are introduced for the new variables by commuting 
space and time derivatives. This analysis has the drawback that the resulting 
system is considerably enlarged and many new constraints appear, although in 
|14) the propagation of these new constraint is trivial. Another approach, the 
one sketched here, is to add as new variables the square root of the Laplacian 
of some of the original variables and so get a first order pseudo-differential 
system. The advantage of this method is that now the principal part of the 
pseudo-differential system is algebraically much simpler to deal with, and the 
new constraints arising from the introduction of the new variables are fewer 
and their propagation is trivial. This procedure has been used by Taylor to 
describe higher order operators, in a general setting and by JI] for Einstein's 
equations as a fully second order system. We shall consider here only constant 
coefficient systems. To deal with constant coefficient systems is all what it is 
actually needed, for if they are well possed, then well possed holds for all systems 
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whose "frozen coefficient" versions coincide with them. 

We consider the ADM equations Unearized off Minkowski space-time, 

dtShij = ~2kij 
dthj = ift.''' [-dkdiSh.j - (1 + h)didj5hki + 2c'fe%(5/ij),] 

where before hnearizing the lapse has been densitized as iV = {g)''Q where Q is 
considered a given function of space-time and the shift has been set to zero. 

We now introduce a new variable, iij = ^/A^Shij, and the corresponding 
evolution equation: 



dAj = -2\/A% (9) 

Substituting in the second equation of jnj all hij for the new variable we get 
a first order pseudo-differential system of evolution equations. The necessary 
and sufficient conditions for well posedness of this type of systems are the same 
as those of differential systems. In fact the proof is the same. Thus it is enough 
to look at the principal symbol, [Namely the limit of the right hand side of the 
equation where all derivatives have been substituted by iLUj divided by |w| when 
\lu\ — > oo, and subsequently multiplied back by |ti;|.] which in this case is: 



A 



% ) \ -(1/2) (4 + (1 + b)Cj,Ujh'''iM - 2Cj''u^ilj)k) 



where we have substituted each partial derivative by cji, and defined, uji := 

UJi/\uj\. 

Thus what we need to know is whether matrix A can be diagonalized. It 
turns out that it can not, although all eigenvalues are real. So the system is 
only weakly hyperbolic. Furthermore it can be shown that the addition of the 
Hamiltonian constraint does not resolve the degeneracy (although it changes 
the characteristic speeds and the subsidiary system of equations for constraint 
propagation). 

The BSSN equations are a further generalization of the above system in 
which some linear combination of first derivatives of the metric are promoted to 
new independent variables and in the corresponding new evolutions equations 
the momentum constraint is conveniently added. The new variable is defined 
with respect to a background connection di and is given by fi :— \hf^^{dkhii — 
dihki) and in its evolution equation the momentum constraint is added (with a 
new parameter c) as follows: 

dth :- -29fefc^ + d,k + c{dkk\ - d,k) (10) 
In this case the principal symbol becomes: 
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A 



( 




) 



\u\ -(1/2) (4 + bCo,u,h>''iki - 2Co(Jj) 



) 



It turns out in this case that for non-zero values of c the symbol is diagonalizable 
with all eigenvalues real, so the system is strongly hyperbolic in the algebraic 
sense. Furthermore one can show that the symmetrizer is a smooth function of 
all variables, so that the well posedness proof for the fully quasi-linear, pseudo- 
differential case applies. This we resume in the following 13 : 

Proposition 6 The BSSN family of systems with c > is well posed. 

Thus we have seen that the application of the concept of strongly hyperbol- 
icity, although used here in the broader class of pseudo-differential quasi-linear 
systems facilitates the well posed analysis of higher order evolution systems. 

5 Constrained Systems 

In this section we deal with another issue for which the notion of strong hyper- 
bolicity is critical. That is the initial value problem for constrained systems. 

In physics one usually finds that the fields under consideration do not only 
satisfy some evolution differential equations but also are subject to some extra 
differential relations, called constraints, which can be thought of as relations at 
a given time, thus effectively diminishing the degrees of freedom of the theory. 
In order for a theory of this type to have solutions some integrability conditions 
between these two set of equations must be satisfied. Part of these integrability 
conditions basically assert that if the constraint equations are satisfied at some 
initial surface, then their time derivative off such surface must also vanish, that 
is, these derivatives should be expressed as homogeneous functions of the con- 
straints themselves and their derivatives tangent to the initial hypersurface. We 
shall call these evolutionary integrability conditions, the subsidiary system; they 
are equations on the constraint quantities, namely the quantities representing 
the constraint equations. If the evolution equations as well as the constraint 
equations are first order in derivatives, something we shall assume from now on, 
then these evolutionary integrability conditions, form also a semi-linear system 
of first order differential equations on the constraint quantities. 

The evolutionary properties of these subsidiary systems are very important. 
At the continuum level in asserting whether or not the constraints will be pre- 
served (continue to hold) along evolution. In particular, when boundaries are 
present, not all otherwise allowed boundary conditions are consistent with con- 
straint evolution. In order to elucidate which ones can be given it is usually 
necessary to have a detailed knowledge of the properties the subsidiary system. 
But even more at the discrete level, for in numerical simulations the constraints 
are never exactly satisfied and so it is necessary to study how these errors will 
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propagate along the simulation. Thus, it is natural to ask which properties, if 
any, these subsidiary systems inherit from the evolution equations. This is going 
to be the main task of this section. We shall see -under certain assumption on 
the constraint system structure- that if the evolution system is strongly hyper- 
bolic then so is the subsidiary system, and furthermore the characteristics of the 
subsidiary system are a subset of the characteristics of the evolution system. 

We next define the set of equations we are going to consider. For a more 
formal definition of many of the concept here introduced see jHI- 

Deflnition 5 By an evolution system with differential constraints we mean: 

i) a first order evolution system of quasi-linear partial differential equations 
of the form, 

dt^ = ^(0, t, xy'^pd.^p + B{^, t, xr, (11) 

ii) a first order system of quasi-linear partial differential equations -the con- 
straint system- of the form, 

:= K{(t>, t, xf^pd^q}" + L{(t>, t, x)"^ = 0, (12) 

Hi) an identity derived from the above two systems which we shall call the 
subsidiary evolution equation. It is a subset of the integrability conditions, and 
we assume has the form, 

dtC^ = S{<f>, t, xy^Bd^C'' + R{c^, d(j), t, x)^bC'', (13) 

Condition Hi) ensures that, if vanishes at the initial surface, hen its 
derivative off the initial surface also vanish. This is clearly a necessary condition 
for the constraints to be satisfied away from the initial surface -it is part of the 
integrability conditions- but in general it is not sufhcient. If the integrability 
conditions exist for a first order quasi-linear systems, then a subset of them 
must have the above form. In particular, if a further condition is assumed (to 
be discussed below) it must be a linear equation on the constraint quantities, 
but that does not play a role in the present discussion. The above identity 
should arise by taking a time derivative of the quantities and using the 
evolution equation for 0". Then re-expressing all quantities in terms of C^. If 
that is so we then say that the system is constrained and are the constraint 
quantities. Failure of the identity implies the existence of further constraints. 

Equation l(T^ is in general not unique, and so S{(j),t,x)^°' b- Thus, in this 
generality, nothing can be said and there are simple examples which show that 
some subsidiary systems do not inherit any good property form the evolution 
systems from which they arise. This happens if there is some differential identity 
between the constraint quantities. That is, if there exists W"" a such that 

iy"AV„C^ = F(C^). (14) 
Since Eq. ({Tl|l is an identity for arbitrary values of 0" we must have that 
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W^^'aK^'^K = 0. (15) 
The following example 4 illustrates the problem: 

Example 2 Consider the evolution equation, 

9*0" = ra,</)" (16) 

in flat space-time where /* is any non-vanishing vector. Let the constraints he: 

C," = (17) 
Then Ki"'"'^ = S'^pS^'i. Since partial derivatives conmute we have 

%c,] = d[3d,^r = (18) 

Thus, any field satisfying W"'^ p = W^"'^'^ p. Gives rise to an identity, and in 

fact all of them satisfy Eq. W^"^ aK^'^K = ^^^"'^'if/^^a = M^a"™^ = 0. 

Thus, the set of subsidiary conditions is given by: 

atQ" = vd,^^ + MA"'^"'0a„c/ (19) 

where W"'"'"'/? is arbitrary. Taking, for instance, VF"'"'"'/3 = we get a sym- 
metric hyperbolic system, hut taking instead, VF"'"'"';3 = —25°'fj5i^l^'^ we get a 
system which is only weakly hyperbolic. 

To exclude this indeterminacy we require the following condition: 
Assumption: For any Wa, and any (j)" , K{(f),xy^ i^uji is surjective. 

This ensures there is no nontrivial W"a satisfying (|15|l . This condition 
ensures that all constraints are differential, as opposite to merely algebraic ones, 
and there are no superfluous ones that can be deduced from any differential 
identity. This property also ensures that the subsidiary system is linear in the 
constraint quantities. 

Algebraic constraints are dealt in a different way beforehand (basically by 
restricting the set of variables to those satisfying them), so in that respect 
this is not a real loss of generality. But there are many important constrained 
systems which do not satisfy this condition, although in many cases one can 
choose a subset of constraints which do satisfy the condition, and the rest has a 
trivial evolution. That is the case in the two examples (representing Maxwell's, 
and Einstein's systems respectively) given below. There are cases in which 
there are sets of measure zero for which the surjectivity property fails, [In the 
example above instead of taking Ci" as the set of constraints, just take the subset 
C" := s'^Ci", for arbitrary constant vector s'.], those have to be treated in grater 
detail, for it could happens that for those values one looses the properties one 
is looking for. 
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Example 3 Maxwell's equations 0) as introduced in the first section, for the 
electric field and vector potential. Two terms proportional to the constraints 
have been added. They are parameterized respectively by the numbers a and [3. 
In order for this system to represent only Maxwell's equations constraints must 
be imposed, indeed the full constraints needed to make the system equivalent to 
Maxwell's equations are given by Cijk = 2d[iWj]k = 0, C — d''Ei — 0. This set 
of constraints does not satisfy the surjectivity condition, but there is a subset 
of them, given by (C := d^Ei,Ci :— d^Wij — diWk'^), whose K^^ satisfies the 
surjectivity condition and furthermore has a closed integrability condition. 

dtC = -ad'C^ 

dtC ^ -PO^C. (20) 

Once this subset is shown to vanish, we can use the subsidiary equation 
dtCijk — Pe^idj]C to ensure the vanishing of the whole set of constraint equa- 
tions. 

Example 4 The following symmetric hyperbolic system describes linearized gen- 
eral relativity in flat space-time ^ , fl^ . JBl 

dthj = -e^'dkfHj 

dtfkij = -dkhj (21) 

where k^j = k^^j-^, and fj^^ij — fk{ij), as before e'^ is the the flat Euclidean 
metric. It admits a one parameter family of constraints: (Ci :~ (kij — 
e.jfcfc'=),Cy d'^f^.k - djfu' + a{dj(fa' - fi\) - e,,{d^{fj - /,'„)). For 
a ^ ^ the corresponding linear map K^^a is surjective, and the subsidiary 
system is given by: 

dtCi = —d^Cki 

dtC,, = -d,Cj+a[d,C,-e,jd''Ck) (22) 

This system has more constraints, but they are preserved (satisfy an homo- 
geneous ordinary evolution equation), if the above ones vanish (see 13}). 

We now focus on the main result of this section. Taking a time derivative 
of (|ll|l and using \i2\ , it is a short calculation to show that for H13|l to be valid 
for arbitrary 0" the following algebraic identity must hold: 

^(<|A|^^j>^ _ Si^\MgK^)B^ ^ (23) 

From this identity we can deduce the following results: 

^The analysis for the non-linear system is identical, since the frozen coefRcient principal 
symbol is the same (up to the addition of trivial shift contribution and multiplication by the 
lapse). 



17 



Lemma 1 Given any fixed non-vanishing co-vector Ua- If{a,(j)°') is an eigenvalue- 
eigenvector pair of A^^f^LOi then (cr, w"^ — K'^^^Wiu"), if is non-vanishing, is 
an eigenvalue-eigenvector pair of S^^b^i- 

Proof 4 Contract the above identity twice with uji and then with . 

Proposition 7 If the above assumption holds, then if A^°' puii has a complete 
set of eigenvectors so does have S^^b^i- Furthermore the set of eigenvalues of 
S^^B^i is a subset of the set of eigenvalues of A^°' pWi. 

Proofs If the eigenvalues of A^°'fjUJi, {u"} expand the whole space, then the 
set {v^ = K^^a'^iu"'} expands the rank of K^'^ai^i, but this map is assumed 
to be surjective, so they expand the whole space. Since they are all possible 
eigenvectors, and since they have the same eigenvalues as the the set {u"} from 
which they are built on, we conclude that all the eigenvalues of S^^B^^i are also 
eigenvalues of A^"' puji. 

In particular we have the foUowing trivial consequences: 

Corollary 1 If A^°' puji has a complete set of eigenvectors and if all their eigen- 
values are real, then so are the eigenvalues of S^^ b^i- 

Corollary 2 The characteristics of the subsidiary system are a subset of the 
characteristics of the evolution system. 

The above Proposition and its first Corollary imply the main result: 

Proposition 8 // the evolution system is strongly hyperbolic, then so is the 
subsidiary system. 

Remarks: 

• By strongly hyperbolic here we refer to the algebraic properties, we have 
not explored the possibility that with these conditions one could conclude 
that if the symmetrizer of the evolution system is smooth then the sym- 
metrizer of the subsidiary system would also be smooth ^ . That might 
very well be the case. 

• There remains the question as to whether a symmetric hyperbolic system 
would yield a symmetric hyperbolic subsidiary system or the converse, 
namely if a non-symmetric but strongly hyperbolic system could give rise 
to a symmetric hyperbolic subsidiary system. 

In example 0| the system is strongly hyperbolic for all values of the pa- 
rameter a. But there are no symmetrizers when a < ^ or a > 1 for the 
subsidiary system. Thus a symmetric hyperbolic system gives a subsidiary 
system which is only strongly hyperbolic. 

^We refer to the first section for further comments on this point. 
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In example 121 the two parameter family system is strongly hyperbolic for 
all values of the parameters such that a(3 > 0. But no symmetrizer exists 
for a > or /3 > Thus they are not symmetric hyperbolic. Never- 
theless, the constraint quantities, have a symmetric hyperbolic subsidiary 
system for all values of the parameters for which the above one is strongly 
hyperbolic. 
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